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INVESTIGATION OF REAL-GAS AND VISCOUS EFFECTS ON THE 


AERODYNAMIC CHARACTERISTICS OF A 40° HALF-CONE 
WITH SUGGESTED CORRELATIONS FOR THE 
SHUTTLE ORBITER 

James L. Hunt, Robert A. Jones, 
and William C. Woods 
Langley Research Center 


SUMMARY 

The purpose of this paper is to evaluate parameters that might be used 
to correlate shuttle orbiter aerodynamic data to be used in extrapolating 
from wind-tunnel to flight conditions. Preliminary calculations indicate 
that the lee-side forces will have an insignificant influence on the aero- 
dynamic characteristics of the orbiter for moderate angle-of-attack entries; 
therefore, this work is focused on phenomena which have an overriding influ- 
ence on windward forces, namely, real-gas (equilibrium and nonequilibrium) 
and viscous-interaction effects. Analytically determined flow fields pre- 
viously obtained on 40^ blunted cones were used as a data source to evaluate 
various correlation parameters. Inviscid effects were found to be the domi- 
nant contributor to the aerodynamic coefficients (40^ half -cone) in the alti- 
tude range of 64 to 76.2 km. The most suitable correlation of the aerodynamic 
forces on these cones is based on local dynamic pressure and local Mach number. 


INTRODUCTION 

As a blunt body at moderate to high angles of attack enters a planetary 
atmosphere at hypersonic speeds, the gas molecules processed by the shock are 
excited to higher vibrational, chemical, and ionization energy modes, and absorb 
energy from the post-shock flow field. As additional energy is absorbed, the 
conservation laws and the thermophysics of the gas dictate certain changes in 
the shock-layer flow (ref. 1). The static temperature, speed of sound, and 
velocity in the real-gas shock layer are reduced; the density is increased and 
the shock-layer thickness is reduced in proportion to this increase provided 
dissociation is not driven near completion. The energy transfer from the shock- 
layer flow to the molecular energy modes and vice-versa downstream of the shock 
occurs at finite rates, which can be characterized by relaxation times, or by 
relaxation lengths in a flow of given velocity. If all relevant relaxation 
lengths are very much shorter than the smallest flow-field dimension of inter- 
est, the flow is regarded as being in thermochemical equilibrium (ref. 2) and 
the greatest departure (largest real-gas effects) from the "ideal” (nonreacting, 
constant ratio of specific heat y) post-shock flow field occurs when energy 
is being transferred to the molecular modes. If all relevant relaxation 
lengths are much greater than the largest flow-field dimension of interest, 
the flow may be regarded as being frozen (ref. 2) so that no energy transfer 


between molecular modes and the post-shock flow field occurs. In this situa- 
tion the flow media may be treated as an ideal gas with constant specific 
heats. If relaxation lengths and flow-field dimensions are comparable, depar- 
ture from thermochemical equilibrium will occur (ref. 2). The resulting non- 
equilibrium (finite-rate) real-gas effects on aerodynamic properties should 
lie between equilibrium and frozen flow provided ionization is not present. 

Hypersonic wind-tunnel aerodynamics data (ideal gas) (pitching -moment 
coefficient Cm, drag coefficient Cj), and lift-drag ratio L/D) on slender 
bodies at small angles of attack have been correlated with some success on 
the basis of a hypersonic flat-plate viscous-interaction correlating 
parameter V* 


V V = 


1 

OO , (• 

where the Chapraan-Rubesin free-stream constant is 

^ CO 

° 


( 1 ) 


and 

M„ free-stream Mach number 

^ 00,1 free-stream Reynolds number based on local conditions 

y’ viscosity evaluated at reference temperature 

free-stream viscosity 
reference temperature 
free-stream temperature 

The correlation of C^, Cj^, and L/D in terms of for orbiter wind-tunnel 

data has been explained on the basis of a change in skin-friction drag. (See 
ref. 3.) The questions are what effect will the real-gas conditions of flight 
have on 0^, Cj), and L/D and to what extent will these effects influence 
the correlation of the aerodynamic data. Along with attempting to extend this 
correlation to include real-gas effects, several other factors should be con- 
sidered. First, since viscous drag is sensitive to local flow conditions at 
the edge of the boundary layer and since these local flow conditions are 
changed in flight because of real-gas effects, local flow conditions should 
be used. Second, the effect of other factors such as lee-side pressure and 
real-gas effects on the windward-side pressure coefficient and thus on 
Cj), and L/D should be considered. 
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The guidelines for the task of devising credible correlating parameters 
are based on ( 1 ) the relative contribution to the vehicle arerodynamics of the 
influencing phenomena (whether inviscid, viscous, or interaction) for a given 
portion of the entry trajectory and (2) the fact that real-gas effects result 
in large changes in local boundary-layer edge conditions; therefore, any attempt 
to correlate aerodynamic data should begin by using local conditions to evaluate 
correlating parameters. 

Viscous interaction parameters are derived by using local conditions at the 
edge of the boundary layer (ref. 4) which account for changes in the edge condi- 
tions from that of the free stream due to the vehicle angle of attack (shock 
processing of the flow) and accompanying inviscid real-gas effects. However, 
evaluating the viscous interaction parameters at the edge of the boundary layer, 
regardless of which weighted reference temperature is used in evaluating the 
Chapman-Rubesin constant C*, does not account for viscous real-gas effects 
(dissociation and recombination in the boundary layer) which may greatly alter 
boundary- layer profiles and thus vehicle skin friction and induced pressure. 

One way to include viscous interaction in the evaluation of orbiter real- 
gas effects on boundary-layer edge condition, as well as those that occur within 
the boundary layer due to dissociation, is to employ a partially coupled 
inviscid-viscous real-gas solution (equilibrium and finite rate) to calculate 
viscous and inviscid components of drag and pitch and correlate the results. 
Since these calculations for the shuttle geometry are a formidable task, the 
approach here was to use a simplified shape (representative of the orbiter con- 
figuration at angle of attack) for which exact solutions can be obtained. This 
has been done on a blunt 40^ half-cone for altitudes from 64.0 to 75.2 km and 
is reported herein. 


SYMBOLS 

A surface area 

Apf> element planform areas (fig. 1) 

^tpf total planform area (fig. 1) 

C’ Chapman-Rubesin constant based on reference temperature, 

y(T’ )T/y(T)T» 

C* Chapman-Rubesin constant based on reference enthalpy, 

M(T’)h/M(T)h’ 

^D,f drag coefficient due to skin friction (viscous drag) 

Cd,p drag coefficient due to pressure 

CD,t total drag coefficient 
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f 

Cm,fD 

^m,pD 

t 

Cp 

D 

h 


h» 


L 


1 

ic 

^ref 

Vpf 

W 

M 

mi 

P 

Pt,2 

Q 

R 

^pf 


pitching-moment coefficient 

moment coefficient due to skin friction 

moment coefficient due to skin-friction drag 

moment coefficient due to pressure 

moment coefficient due to pressure drag 

total pitching-moment coefficient 

pressure coefficient 

drag 

enthalpy 

reference enthalpy, 0.5h^ + 0.22haT^ + 0.28h 
lift 

length from nose apex to local position 
cone axial length (fig. 13) 
reference cone length (table III) 
reference planform length (fig. 1) 
vehicle length 
Mach number 

element moment arm to center of gravity (fig. 1) 
static pressure 

total pressure behind normal shock 

dynamic pressure 

Reynolds number 

cone base radius 

cone nose radius 

reference planform area (fig. 1) 
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^ref 

T 


Tf 


u 

V’ 


X 

y 

z 

a 

Y 


one-half cone base area 
static temperature 


reference temperature, 


T 


1 + 0,0225 



+ 


^0.695 



velocity 

viscous interaction parameter, 
axial distance 

distance normal to cone surface 
positive normal-force direction 
angle of attack 
ratio of specific heats 



d log p 

^ isentropic exponent, 

d log p 

6 flow deflection angle 

0 cone half -angle 

y viscosity 

p density 

(j) cone-ray meridian angle 

Subscripts: 

aw adiabatic wall 

Eq equilibrium 

e boundary-layer edge 

FR finite rate 

1 based on local conditions (at x/i^ = 0.5 for cone analysis 

(fig. 13)) 



Iq based on eone length 

coefficient (drag or pitch) nondimensionalized by local 
dynamic pressure at x/Xq = 0.5 

t total 

w wall 

2 based on local conditions behind oblique shock 

“> based on free-stream conditions 

Abbreviations: 

CFHT continuous flow hypersonic tunnel 

Eq equilibrium 

FR finite rate 

VDT variable density tunnel 

Primes denote quantity evaluated at reference temperature. 


ROUGH ESTIMATES OF EFFECTS OF FLIGHT CONDITIONS 
Leeward Pressure Effects 

To obtain a rough order of magnitude estimate of the effect of lee-side 
pressure on orbiter pitching moment, an estimate was made by assuming con- 
stant pressures to act over the plan view with the geometry of figure 1 . 

With zero pressure coefficient for the lee side at an angle of attack of i|0° 
and Cp = 1.05 on the windward side (y = 1.^ ideal-gas oblique shock value 
for M^ = 20), the pitching-moment coefficient C„j is -0.54. Changing the 
lee side Cp to ±0.03 rather than zero gives only a 3-percent change in C^. 
Therefore, it is felt that any M^, R_^, and real-gas effects would have only 

a small effect on the lee-side contribution to Cj„. 

Furthermore, one estimate of real-gas effects on lee-side pressure indi- 
cates that ideal gas air data may be applicable to flight provided the correct 
Mach number and Reynolds number are simulated. This estimate was on the basis 
of finite-rate chemistry calculations made under contract for shuttle condi- 
tions (ref. 5). They were given an assumed pressure distance variation for 
expansion from the windward to the leeward side (fig. 2) and a flight condition 
of 6.1 km/sec at an altitude of 67.1 km, and a n 20°. The calculated results 
for both equilibrium and nonequilibrium (finite-rate) real air are shown in 
figure 3. On the windward side, at the start of the expansion the isentropic 


6 



exponent is 1.12. If the real gas were to expand around the vehicle and remain 
in equilibrium, the isentropic gamma (d log p/d log p) remains constant at 1.12 
as indicated by the slope of the curve in figure 3 ; however, the finite-rate 
chemistry calculations showed that for the data of figure 2 , the gas would 
actually "freeze” and follow the nonequilibriura (finite-rate) line shown in fig- 
ure 3 which has a slope, and thus a gamma, of 1.43. This result indicates that 
to best match the pressure/density relation in an expansion to the lee side in 
flight at this one condition, one should test in an ideal gas having a gamma 
of 1.43 at the correct Mach number and Reynolds number. In other words, con- 
ventional air wind tunnels may give the best simulation for this one flight 
condition. 


Real-Gas Effect on Windward Pressure Coefficient 

A rough estimate of the effects on due to changes in Cp on the 

windward side of the orbiter was made by using the constant-pressure planform 
approach shown in figure 1 . The magnitude of real-gas effects on Cp for 
oblique shocks is shown in figure 4. Here Cp is shown for several ideal 
gas values of gamma at Mach 20 (oblique shock calculations for gamma of 1.1, 

1.4, and 1.6?) and for an equilibrium real gas (realistic for orbiter windward 
wide) at an altitude of 60 964 m and a velocity of 6096 m/sec (ref. 6 ). For 
small deflection angles real-gas effects are small, but they become large as 
the flow deflection angle increases. For an oblique shock at 6 = 40^^ which 
corresponds to a = 40^, the real-gas Cp is 0.95 whereas the ideal-gas air 
value at M^ =20 is 1.05. Integrating as shown in figure 1 indicates that 
the ideal gas Cjjj is -0.054 whereas the equilibrium real-gas Cj^ is -0.049, 
a 10 -percent change. 

A better estimate of the real-gas effect could be made by using wind-tunnel 
test data at various effective gammas. The effective value of gamma (isentropic 
exponent (ref. 6 )) in the shock layer for flight varies over the trajectory as 
shown in figure 5. This value varies from about 1.1 to 1.24 during the high- 
velocity high a part of entry. The Langley hypersonic CF 4 tunnel (ref. 7) 
is capable of obtaining force and moment data on a model of the orbiter at 
these effective gammas. 

Questions arise as to how gamma varies over the windward side in the shock 
layer and over the trajectory. Figure 6 may help to show the relative insen- 
sitivity of gamma for the orbiter flight conditions. This figure is taken from 
reference 8. The value of y remains approximately constant for a wide range 
of velocities and is only weakly affected by pressure for orbiter hypersonic 
entry conditions (matches CF 4 tunnel). 


Real-Gas Effect on Elevon Effectiveness 

The downward deflection of the elevon of the orbiter during hypersonic 
flight will generate an embedded shock. The inviscid real-gas effects on Cp 
ratio for this condition are shown in figures 7 (flight fairing from ref. 9 ) 
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and 8 and compared with that for ideal-gas wind tunnels. As shown, the 
presence of the embedded shock causes large differences in the eleven pres- 
sure ratio between real-gas and wind-tunnel conditions. 

Another factor which complicates this analysis is the possibility of 
separation of the low Reynolds number boundary-layer flow .due to the adverse 
pressure gradient. Extensive boundary-layer separation could result in greatly 
reduced elevon effectiveness. 


CORRELATION OF ORBITER HYPERSONIC AERODYNAMIC DATA 
Similitude Parameters 


pjeal-gas effects result in large changes in local Mach number behind 
oblique shock M 2 , local Reynolds number R 2 , as well as in C*. The resulting 
changes in the correlating parameters 


V2 


M2\|^ 



( 2 ) 


are shown in figures 9 to 12 and compared with the values for an orbiter tra- 
jectory as well as the values that can be obtained in several wind-tunnel facil- 
ities . Obviously, for either or V 2 , if Mach number and Reynolds numbers 

are matched, V’ is simulated rather closely. However, just matching V’ does 
not mean that Mach number and/or Reynolds number are matched. As can be seen 
by comparing figures 9 and 10, it is possible to test at flight or V 2 

and be in error in both Mach number and Reynolds number. Since aerodynamic 
similitude in ideal gas requires the duplication of M^, R^, and T^/T|-, the 

question becomes that of the degree of coupling between these so-called inde- 
pendent parameters. Also, the question as to which forces are dominant (invis- 
cid, viscous, or interaction) must be addressed before a correlation parameter 
can be selected for the total (inviscid or viscous) aerodynamic characteristics 
of a vehicle. 


The viscous interaction parameter may be a good parameter for corre- 

lating ideal- gas hypersonic aerodynamic data for a given y test media; how- 
ever, it is probably a poor correlation parameter for flight-entry hypersonic 
aerodynamic data for the following reasons: (1) Viscous effects are dependent 

on the boundary-layer profiles which are, in turn, influenced by the chemical 
reaction in the boundary layer. A viscous interaction parameter cannot account 
for these influences, (2) Viscous effects are governed by the local conditions 
at the edge of the boundary layer (not free stream) and real-gas effects greatly 
alter these local conditions in flight as compared with those in ground facil- 
ities (figs. 11 and 12). A viscous interaction parameter based on free-stream 
conditions cannot account for these alterations. (3) The real-gas influence 
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on pressure level and distribution probably has a larger effect on the aero- 
dynamic coefficients than viscous interaction over a significant part of the 
entry trajectory. 


Flow-Field Solution Technique 

Since aerodynamic center of pressure shifts due to real-gas effects and 
viscous interactions along with changes in control surface effectiveness may 
alter the trim condition of the vehicle, it is imperative that the shuttle 
orbiter aerodynamic correlation parameters be properly evaluated. One means 
of evaluating these real-gas and viscous interaction effects is to employ a 
partially coupled inviscid-viscous real-gas solution to calculate viscous and 
inviscid components of drag and pitch on the orbiter, but this is a formidable 
task; therefore, one approach is to make such calculations on a simplified con- 
figuration at orbiter flight conditions and correlate the results. Such cal- 
culations have been made on a blunt 40^ cone. The flow-field calculations are 
axisymmetric but the aerodynamic coefficients were obtained by integrating over 
one-half the cone. Viscous and inviscid components of drag and pitch were 
obtained by using quasi-coupled inviscid-viscous real-gas solutions of a 
40^ half-angle cone (fig. 13) for both orbiter flight conditions (altitudes 
of 64.0 and 76.2 km) and representative tunnel tests. 

The flow-field solutions (equilibrium and finite-rate) required in obtain- 
ing the flow-field results presented in tables I to IV were calculated by using 
previously developed computer codes to define inviscid and viscous flow fields. 
The computer codes utilized (NASA Ames equilibrium blunt body and method of 
characteristics, modified Curtis and Strom unified nonequilibrium flow field, 
and a modified version of the viscous reacting gas code developed by Blottner) 
and the manner in which they were applied are given in reference 10. The pro- 
cedure utilized in reference 10 was not fully coupled in that the inviscid phase 
(inviscid flow-field properties) and viscous phase (boundary layer) were com- 
puted separately; however, these two independent phases were interfaced through 
a viscous-inviscid mass flux matching in the shock layer along points at the 
edge of the boundary layer. Initial estimates of the edge conditions for the 
boundary layer were made and the computation carried out. The resulting 
boundary-layer mass flow was then used to interpolate along the previously 
computed inviscid field rays to determine flow conditions in the inviscid field 
at a point corresponding to the boundary-layer mass flow value. The results 
herein are thus subject to the qualification that the boundary-layer displace- 
ment thickness effects (displacement thickness corrections to the cone surface 
contour) are small, and at least for the cases analyzed (tables I to V), the 
displacement thickness correction need not be considered in the inviscid flow- 
field computation. 


Objective 

The objective of this effort is to correlate viscous interaction and real- 
gas effects for orbiter-like entry and tunnel conditions (analytical, no experi- 
ment) for the purpose of extrapolating tunnel aerodynamic data to flight condi- 
tions. The assumption herein is that parameters which correlate the aerodynamic 
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data on this simplified configuration (40° half -angle cone) should also be 
applicable to the shuttle orbiter at entry conditions; therefore, the emphasis 
of this paper is on selecting an appropriate set of correlation parameters for 
the 40° half -angle cone aerodynamic coefficients. (See tables I to V.) Such 
correlations are critical to the successful prediction of the orbiter aerody- 
namic characteristics during entry since aerodynamic center shifts due to real- 
gas effects not yet accounted for and viscous interaction improperly correlated 
along with changes in control surface effectiveness may alter the trim condi- 
tions of the vehicle. 

The guidelines for this task of devising correlating parameters are based 
on (1) the relative contribution to the vehicle aerodynamics of the influencing 
phenomena (whether inviscid, viscous, or interaction) for a given section of 
the entry trajectory and (2) the fact that real-gas effects result in large 
changes in local boundary-layer edge conditions. Therefore, local conditions 
must be considered in evaluating correlating parameters. 


RESULTS AND DISCUSSION 

Both inviscid and laminar viscous (coupled) real-gas (equilibrium 
and nonequilibrium) flows were calculated for a series of blunt cones 
(0 = 30° and 40°, orbiter flight conditions 

(64.0 km ^ altitude = 76.2 km, 4-9 = velocity = 7.3 km/sec) as well as for 
ideal-gas conditions in several wind tunnels (ref. 10). The skin friction 
and pressure were integrated over one-half of the surface of the 40° half- 
angle cone (ic/i^n = 13) for the various cases already computed 

(ref. 10) along with computing several local parameters for correlative 
purposes. The results are given in tables I to V. 


Viscous Interaction Parameters 

Since hypersonic wind-tunnel aerodynamic data on slender bodies at small 
angles of attack have been correlated with some succes^ on the basis of the 
hypersonic viscous interaction correlating parameter V^, the viscous coeffi- 
cients on the 40° half-angle cone (fig. 13) will be first examined on the 
basis of this parameter. Figures 14 and 15 show the viscous drag and 
pitching-moment' coefficients as a function of for both the flight and 

tunnel calculated data on the 40° half -angle cone. The data correlate well 
on a straight line except for the Mach 20 nitrogen tunnel data_. The same 
coefficients are shown in figures 16 and 17 as a function of V based on 
local coefficients at x/1q = 0.5. The data for both viscous drag and 
pitching-moment coefficients including the Mach 20 nitrogen tunnel cases 
correlate along a straight line and indicate What a better correlation of 
the viscous coefficients is obtained with a V based on local conditions 
(x/Iq = 0.5) rather than on free-stream values. 

If local conditions are to be used in the correlation, then the aerody- 
namic coefficients themselves should be formulated with a local dynamic pres- 
sure rather than with the free-stream dynamic pressure. Thus, the viscous 
drag coefficient based on the local dynamic pressure external to the boundary 
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layer at an x/Iq =0.5 is given in figure 18 as a function of V^. Here 
the flight and CF4 tunnel data (ref. 10) correlate along one straight line 
whereas the remaining wind-tunnel data correlate along a different straight 
line. There are probably two reasons for this difference: (1) local Mg, Rg, 

and T^/Tt are independent simulation parameters and they are not coupled into 
one independent simulation parameter by and (2) boundary-layer edge con- 

ditions are similar in flight and in the CF4 tunnel. The second stems from the 
fact that viscous effects not only depend on the boundary-layer edge conditions 
but on boundary-layer profiles which are in turn influenced by the chemical 
reaction in the boundary layer. This result is supported by the boundary-layer 
temperature profiles shown in figure 19. The profile shown for CF4 at = 6 
is much more similar to the two flight profiles shown than those calculated in 
air, helium, and nitrogen hypersonic tunnels. 


Relative Contribution of Phenomena Influencing 
Aerodynamic Coefficients 

Before attempting a correlation of the total (inviscid plus viscous) aero- 
dynamic coefficients, the question that must be considered is what percentage 
of the aerodynamic coefficient is contributed by specific influencing phenomena 
(whether inviscid, viscous, or interaction). The percent of the drag coeffi- 
cient contributed by the viscous, finite-rate, and real-gas effects is shown in 
figure 20 as a function of V^. For flight conditions from an altitude of 64.0 
to 76.2 km, viscous effects account for only 1 to 2 percent of the total drag 
coefficient. Finite-rate effects are even smaller. The differences in tunnel 
and flight drag coefficients (fig. 20) indicate that the real-gas effects in 
flight account for approximately 10 percent of the total drag. 

The percentage contribution to the pitching-moment coefficient is given 
in figure 21. The viscous effects account for 4 to 10 percent of the total 
pitching-moment coefficient, finite-rate effects again being smaller. Real-gas 
effects are on the order of 5 to 1 1 percent or essentially the same as the 
viscous effects. The obvious implication is that the influencing phenomena of 
the total aerodynamic coefficients for this cone angle in this flight regime 
are 90 to 98 percent inviscid; therefore, the proper correlation must rely 
heavily on inviscid parameters. At higher altitudes or smaller cone angles 
viscous effects could become more significant. 


Data Spread 

The inviscid drag coefficients for both flight and wind-tunnel conditions 
are given in figure 22 as a function of V^. The flight data lie in the lower 
left-hand corner whereas the wind-tunnel data are spread out along the top of 
the plot. Obviously, does not correlate the inviscid drag; however, the 

spread is only approximately 8 percent (real-gas effects). Thus, in searching 
for a more suitable correlation parameter which according to the previous sec- 
tion must rely heavily on inviscid parameters, a fine tuning is required since 
small percentage scatter is involved. Therefore, the criterion for a correlation 
parameter is that the scatter of the data is less than 10 percent. 
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Free-Stream Inviscid Correlation Parameter 

A correlation of the pressure coefficient based on free-stream y and 
M as suggested by the hypersonic similarity law (ref- 11) is given in fig- 
ure 23 . This correlation is not acceptable since the data scatter is much 
greater than 10 percent in the flight part of the curve. Also, the fairing is 
strongly nonlinear which makes extrapolation from wind tunnel to flight condi- 
tions inaccurate. This is very significant here since in an actual extrapola- 
tion of shuttle orbiter wind-tunnel data, flight data points are not available. 
Thus, depending on how the tunnel data are faired (M^ ^ 20), extrapolation 
errors of 20 to 30 percent are possible for flight Mach numbers in the Mach 25 
to 30 range. (See fig. 23. Use of Mach 20 helium data points would make 
tunnel-data extrapolation error much larger.) 


Local Inviscid Correlation Parameter 


Inviscid similitude at geometrically similar points in model and proto- 
type compressible fluid flow systems requires duplication of Mach number and the 
isentropic exponent (gamma); therefore, an inviscid correlation parameter should 
contain both Mach number and gamma. Prior discussion of real-gas effects herein 
also indicates that the correlation parameter should be based on some local ref- 
erence conditions behind the bow shock and not on free-stream conditions. 

It was surmised in the section ^Viscous Interaction Parameters” that if 
local conditions are to be used in the correlation, then the aerodynamic coef- 
ficients themselves should be formulated with a local dynamic pressure rather 
than with the free-stream value. 


Examining the inviscid drag coefficient 


(CD,p)q^ 



I 


p dA 


1 

2 





(3) 


where y 

is now the 

isentropic 

exponent. For a sharp cone with no viscous 

effects, 

(CD,p)q^ 

becomes 




Px \ dA 

2 


(CD,p)q^ = 


(4) 



'h p,M,2a 

I-* 
I— » 


Even though the pressure canceled out (eq. (4)), which it will not do for the 
real situation since bluntness and induced pressure keep the cone pressure p 
inside the integral (eq. (3)), the parameter (C])^p)q^ or (Cm^p)q^ should, 
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when plotted against local Mach number, collapse the gamma effects, 
lation is in essence 


The corre- 



PlA 


YxM^2f(M^) 


(5) 


The pressure drag coefficient based on the local dynamic pressure at 
x/ 1 q = 0.5 is given as a function of local Mach number in figure 24. This 
correlation shows promise since the scatter of the data is less than 3 per- 
cent about a linear fairing of tunnel and flight data points (the Mach 20 
helium data point which is farthest from the flight part of the fairing being 
neglected). Also, the tunnel data fairing is very near linear in the flight 
part of the fairing where the CF 4 tunnel data point lies; this linearity would 
probably allow a more accurate extrapolation of wind-tunnel data to flight 
conditions. 


The inviscid pitching moment based on local dynamic pressure is given as 
a function of local Mach number in figure 25. The data spread about the lin- 
ear fairing of flight and tunnel data is approximately 4 percent and the fair- 
ing of tunnel data is again near linear in the flight part. 

The total drag and pitching-moment coefficients based on local dynamic 
pressure are given as a function of local Mach number in figures 26 and 27, 
respectively. The data spread about the linear fairing in the drag coeffi- 
cient correlation is less than 3 percent. The spread is slightly larger in 
the pitching-moment correlation; however, the correlation is extremely good in 
the flight data region. Also, finite-rate flight data correlated equally as 
well as equilibrium flight data. 


Space Shuttle Orbiter Correlation and Extrapolation 
to Flight Condition 

The correlation of total pitch and drag nondimensionalized by local dynamic 
pressure as a function of local Mach number was extremely good for the case of 
the 40^ half-angle cone in that the data spread about data fairing was less than 
or equal to 3 percent; the tunnel-data fairings were near linear and the CF 4 tun- 
nel data point for both pitch and drag lay within the flight local Mach number 
regime which makes the extrapolation to flight conditions amenable. Also 
finite-rate flight data correlated equally as well as equilibrium flight data. 
This result leads to the question of how such a correlation based on local con- 
ditions is established for the shuttle. 

First, aerodynamic wind-tunnel data are obtained on the orbiter configu- 
ration for a given angle of attack at hypersonic Mach numbers in helium, air, 
nitrogen, tetrafluoromethane (CF 4 ), and hexafluoromethane (C 2 F 5 postulated 
tunnel). Hexafluoromethane has been included because, according to figures 10 
and 12, a local Mach number significantly higher than that which occurs at 
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Mach 6 in CF 4 could be obtained. This condition would result in two experi- 
mental wind-tunnel data points (CF 4 and C 2 F 5 ) in the flight section of the 
pitch and drag fairing (figs. 24 to 27) and thus would lend much more credit- 
ability and accuracy in extrapolation of the correlation. In fact, for the 
altitude range of this study (64 km to 762 km), the addition of a data point 
in C 2 F 5 would alleviate the need for practically any extrapolations. 

Second, the local properties on a cone or wedge with a cone or wedge angle 
compatible with the particular orbiter angle of attack would be calculated for 
each tunnel condition. These, calculated local properties (dynamic pressure and 
Mach number) at an appropriate reference station would be used in correlating 
the orbiter aerodynamic forces. The idea here being that the change in wind- 
ward free stream to local conditions for the orbiter at some reference station 
downstream of the influence of the apex flow could be represented by the change 
in free stream to local conditions of a cone or wedge (with a half-angle equal 
to the orbiter angle of attack) at the same tunnel conditions. A test program 
to establish and validate this correlation for the orbiter in the manner 
described is being pursued. 

Once the correlation has been established with wind-tunnel aerodynamic 
data on the orbiter by using the local Mach number and dynamic pressure calcu- 
lated at the reference location for an appropriate cone or wedge, the local 
conditions on an appropriate cone or wedge at the assigned reference location 
are calculated along the flight trajectory (equilibrium and finite rate). 

These local conditions (equilibrium and finite rate) are used to enter the 
correlation established with wind-tunnel data to obtain the orbiter flight 
aerodynamics. This procedure must be repeated for each angle of attack. 


Aerodynamic Coefficient Correlation for Altitudes Above 
76.2 km and Angle of Attack Lower Than 4QO 

Real-gas effects are extremely sensitive to the angle of attack of the 
vehicle. Decreasing the cone angle below 40° (representing a reduction in the 
orbiter angle of attack) for the 64- to 76.2-km altitude range of this study 
will both decrease the real-gas effects and increase the viscous influence on 
the aerodynamic characteristics. Thus, for some lower angle of attack the 
correlation advocated in the previous section may break down and a V' based 
on local conditions may serve to correlate the aerodynamic coefficients. _If 
the angle of attack and velocity are small enough, a correlation with a V’ 
based on free-stream conditions may be needed. 

At altitudes above 76.2 km for a shuttle orbiter type entry, both finite- 
rate and viscous effects will become more dominant. Here again, the local 
Mach number correlation may collapse but the importance of using local condi- 
tions in viscous correlation parameters will probably become more pronounced. 
One possible solution to the correlation problem as the viscous effects become 
more pronounced is to establish local inviscid correlations such as suggested 
herein for constant levels of viscous effects as indicated by either V{ or 
R,. 
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CONCLUSIONS 


These conclusions apply to a 40^ half-angle blunt half-cone entry along 
an orbiter trajectory in the altitude range from 64 to 76.2 km. Trends indi- 
cated here should be applicable to the shuttle orbiter at an angle of attack 
of 4QO. 

(1) Aerodynamic coefficient (drag and pitch) data spread (tunnel and flight) 
is reduced from an 8- to 10-percent spread to a 3- to 4-percent spread by using 
an inviscid correlation based on local conditions. 

(2) Viscous effects account for 1 to 2 percent of total drag and 4 to 
10 percent of total pitching moment. 

(3) Finite-rate effects account for 0.1 to 2 percent of total drag 
coefficient and 0.1 to 5 percent of total pitching-moment coefficient, 

(4) Real-gas effects are on the order of 10 percent for both drag and 
pitching moment. 

(5) The viscous interaction parameter V’ appears to be a poor correlation 
parameter for hypersonic entry of large angle cones. There are three causes 
for this condition: 

(a) Viscous effects are governed by the local conditions at the 
edge of the boundary layer (not free stream) and real-gas effects 
greatly alter these local conditions in flight as compared with those 
in ground facilities. 

(b) The viscous effects are dependent on the boundary-layer 
profiles which are, in turn, influenced by the chemical reaction in 
the boundary layer. 

(c) A viscous interaction parameter based on free-stream 
conditions cannot account for changes in pressure level and dis- 
tribution which occur because of real-'gas effects. This change in 
pressure has a larger effect on the aerodynamic coefficients than 
viscous interaction over a significant part of the trajectory. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
February 9, 1977 
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TABLE I.- FLIGHT CONDITIONS 


Altitude , 
km 

Velocity, 
km/ sec 

^c/*".l 

0, 

deg 

State 


Free-stream 

conditions 

Local conditions 

(xAo 

□ 0 . 50 ) 

M 

00 

R 1 

OO) Lq 

C 

OO 

OO 





76.2 

7.315 

50 

40 

Equilibrium 

26.1 

0.9540 X 106 

0.4773 

0.01911 

4.038 

0.4050 X 10^ 

1.167 

0.00687 

70.104 

7.315 

50 

40 

Equilibrium 

24.64 

2.0940 

.5055 

.01240 

3.955 

.9350 

1.171 

.00443 

64.008 

6.096 

50 

40 

Equilibrium 

19.50 

3.5014 

.5239 

.00766 

3.762 

1 .8584 

1.162 

.00298 

64.008 

4.877 

50 

40 

, Equilibrium 

15.60 

2.8102 

.5683 

.00713 

3.630 

1 .8141 

1.141 

.00288 

76.2 

7.315 

50 

40 

Nonequilibriura 

26.1 

.9540 

.4773 

.01911 

3.565 

.2600 

1.186 

.00762 

70.104 

7.315 

50 

40 

Nonequilibrium 

24.64 

2.0940 

.5055 

.01240 

3.174 

.4919 

1.177 

.00492 

64.008 

6.096 

50 

40 

Nonequilibrium 

19.50 

3.5014 

.5239 

.00766 

3.762 

1.8584 

1.162 

.00298 

64.008 

4.877 

50 

40 

Nonequilibrium | 

15.60 

2.8102 

1 

.5683 

.00713 

3.630 

1 .8141 

1.141 

.00288 


Altitude , 
km 

— 


0, 

deg 


Reference conditions 

Pressure 

data 


Skin-friction data 

Velocity, 

km/sec 

^c/^n 

State 

kN/m^ 

Vef » 
m 

Sref, 

m^ 

Drag force, 
MN 

Cd,p 

Moment , 
MN-m 

Cm,p 

Drag force, 
kN 

Co,f 

Moment , 
kN-m 


76.2 

7.315 

50 

40 

Equilibrium 

1 .068 

30.48 

1050.55 

0.949 

0.84623 

- 2.991 

- 0.08747 

18.69 

0.01666 

-307.64 

- 0.00900 

70.104 

7.315 

50 

40 

Equilibrium 

2.639 

30.48 

1050.55 

2.348 

.84714 

-7.400 

-.08758 

29.33 

.01058 

-482.86 

-.00571 

64.008 

6.096 

50 

40 

Equilibrium 

4.067 

30.48 

1050.55 

3.634 

.85036 

-11.453 

-.08793 

28.13 

.00658 

-463.13 

-.00356 

64.008 

4.877 

50 

40; 

Equilibrium , 

2.603 

30.48 

1050.55 

2.335 

.85395 

- 7.374 

-.08847 

16.23 

.00593 

-267.20 

-.00321 

76.2 

7.315 

50 ; 

40 

Nonequilibrium 

1 .068 

30.48 

1050.55 

.936 

.83531 

-2.859 

-.08368 

16.36 

.01459 

-269.29 

-.00788 

70.104 

7.315 

50 

40 

Nonequilibrium 

2.639 

30.48 

1050.55 

2.330 

.84058 

- 7.339 

-.08686 

24.38 

.00800 

-401.34 

-.00475 

64.008 

6.096 

50 

40 

Nonequilibrium 

4.067 

30.48 

1050.55 

3.633 

.85027 

- 11.458 

-.08798 

26.64 

.00623 

- 437.03 

-.00337 

64.008 

4.877 

50 

40 

Nonequilibrium 

2. 603 

30.48 

1050.55 

2.335 

.85395 

- 7.374 

-.08846 

15.76 

.00576] 

- 269.29 

-.00311 
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TABLE II.- WIND-TUNNEL CONDITIONS 


Test gas 


1 


Free-stream 

conditions 

Local conditions 

( x/^c 

0 . 5 ) 

00 

^ e , oo , 1(3 

c ’ 

00 



Ri 

D 

n 

Air 

50 

40 

8 

1.2955 X 10 ^ 

0.95172 

0.00686 

2.066 

0.5656 X 10 ^ 

1 .0433 

0.00285 

Air 

50 

40 

8 

7.7728 

.95172 

.00230 

2.066 

3.3920 

1 .0433 

.00114 

Nitrogen 

50 

40 

20 

.0576 

1 .2424 

.09237 

2.208 

.0102 

1.1875 

.02380 

Nitrogen 

50 

40 

20 

.2231 

1 .2424 

.04720 

2.207 

.0396 

1.1875 

.01208 

Helium 

50 

40 

20 

1.9071 

.24187 

.00712 

1 .585 

.1681 

.9800 

.00388 

Freon 

50 

40 

6.2 

.0496 

.95511 

.02721 

3.166 

.0930 

1.0339 

.01056 

Freon 

50 

40 

6.2 

.2978 

.95611 

.01111 

3.166 

.5580 

1.0339 

.00434 


Test gas 

^c/^n 

e, 

deg 

Reference conditions 

Pressure data 

Skin-friction data 

9co> 

kN/m2 

Vef , 

m^ 

Sref , 

m2 

Drag force, 
kN 

Cd,p 

Moment , 
N-m 

^m,p 

Drag force, 
N 


Moment , 
N-m 

B 

Air 

50 

40 

16.591 

0.2134 

0.0513 

0.731 

0.90015 

- 15.519 

- 0.08952 

4.048 

0.00499 

- 0.461 

- 0.00269 

Air 

50 

40 

94.918 

.2134 

.0513 

4.388 

.90044 

- 96.046 

-.09237 

10.186 

.00209 

- 1.166 


Nitrogen 

50 

40 

.807 

.1524 

.0262 

.019 

.90125 

-.298 

-.09077 

1 .201 

.05613 

-.095 

-.03031 

Nitrogen 

50 

40 

3.229 

.1524 

.0262 

.074 

.87397 

1 .101 

-.08784 

2.358 

.02765 

-.190 

-.01493 

Helium 

50 

40 

10.991 

.2134 

.0513 

: .512 

.90744 

10.741 

-.09090 

3.292 


-.380 


Freon 

50 

40 

1.055 

.2134 

.0513 

,049 

.89808 

1 .058 

-.09155 

r .246 

.02272 

-.136 

-.01227 

Freon 

- 50 

40 

6.328 

.2134 

.0513 

.292 

.89818 

6.305 

-.09099 

3.025 

.00093 

-.353 

-.00499 


















TABLE III,- DRAG FORCE MOMENT FOR FLIGHT CONDITIONS 



Reference conditions 

Pressure data 

Skin-friction data 


^ref , 



Drag force 


Drag force 



Sref > 

State 

moment , 

*^m,pD 

moment , 

^m,fD 

kN/m2 

m 



MN-m 

kN-m 


1 .068 

30.48 

1055.55 

Equilibrium 

-10.403 

-0.30422 

185.218 

-0.00542 

2.639 

30.48 

1055.55 

Equilibrium 

-25.733 

-.30456 

290.728 

-.00344 

4.067 

30.48 

1055.55 

Equilibrium 

-26.261 

-.30573 

278.867 

-.00214 

2.603 

30.48 

1055.55 

Equilibrium 

-25.598 

-.30708 

162.454 

-.00195 

1 .068 

30.48 

1055.55 

Nonequilibrium 

-10.223 

-.29920 

162.197 

-.00475 

2.639 

30.48 

1055.55 

Nonequilibrium 

-25.533 

-.30219 

246.461 

-.00292 

4.067 

30.48 

1055.55 

Nonequilibrium 

-26.259 

-.30571 

265.727 

-.00204 

2,603 

30.48 

1055.55 

Nonequilibrium 

-25.596 

-.30708 

157.722 

-.00189 


TABLE IV.- DRAG FORCE MOMENT FOR WIND-TUNNEL CONDITIONS 


Reference 

conditions 

Pressure data 

Skin-friction data 

Q , 


^ref ’ 
m^ 


Drag force 


Drag force 


kN/m2 

^ref > 
m 

Test gas 

moment , 
N-m 

^m,pD 

moment , 
N-m 

Cm.fD 

16.591 

0.2134 

0.0513 

Air 

-55.751 

-0.32181 

- 0.271 

-0.00160 

94.918 

.2134 

.0513 

Air 

-335.999 

-.32314 

-.691 

-.00066 

.807 

.2134 

.0513 

Nitrogen 

-1 .044 

-.32274 

-.054 

-.01815 

3.229 

.2134 

.0513 

Nitrogen 

-4.027 

-.31289 

-.108 

-.00891 

10.991 

,2134 

.0513 

Helium 

-39.102 

-.32475 

-.230 

-.00185 

1.055 

.2134 

.0513 

Freon 

-3.715 

-.32205 

-.081 

-.00735 

6.328 

,2134 

.0513 

Freon 

-22.317 

-.32185 

-.203 

-.00299 
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TABLE V.- ADDITIONAL EDGE PROPERTIES AT x/Xq = 0.5 


(a) Wind tunnel 



Test gas 

Pu 

Ql y 


Ti, 

30 2 

kN/m 

kN/m2 

kN/m2 

L 

m/sec 

K 

16.591 

Air 

14.240 

42.589 

866.55 

437.56 

94.918 

Air 

85.466 

255.345 

866 . 55 

437.56 

.807 

Nitrogen 

.728 

2.480 

1266.75 

819.44 

3.229 

Nitrogen 

2.820 

9.619 

1266.75 

819.44 

10.991 

Helium 

9.973 

20.885 

1671 .22 

320.33 

1 .055 

Freon 

.948 

5.315 

680.01 

436.83 

6.328 

Freon 

5.683 

31.903 

680.01 

436.83 


(b) Flight 


^00 > 


Pi. 

1 

9l. 

Ui , 

Ti, 

kN/m2 

State 

kN/m2 

kN/m^ 

m/sec 

K 

1 .068 

Equilibrium 

0.905 

7.723 

5439.46 

4651 .67 

2.639 

Equilibrium 

2.236 

18.362 

5435.19 

4833.89 

4.067 

Equilibrium 

3.457 

25.439 

4511.95 

3980.56 

2.603 

Equilibrium 

2.222 

15.058 

3598.47 

2986.11 

1.068 

Nonequilibrium 

.900 

6.090 

5407.46 

6162.22 

2.639 

Nonequilibrium 

2.207 

11.960 

5415.99 

5770.00 

4.067 

Nonequilibrium 

3.457 

25.439 

4511 .95 

3980.56 

2.603 

Nonequilibrium 

2.222 

15.058 

3598.47 

2986.11 





^pfl " 

140.25 

mj = 8.5 

%f2 = 

114.4 

mg = -5.2 

%f3 " 

19.24 

mg = -11.7 

A = 

pf4 

2 X 53.04 

m^ = -6.93 

A .c = 

pf5 

2 X 8.17 

mg = -10.97 

^pf6 " 

4.8 

mg = 22.4 

%f7 = 

2 X 5.95 

my = -9.23 


m = Moment arm to center of gravity 
Aj.pj = Total planform area = 396.31 
Spj = Reference area = 271.4 
^rpf “ Reference length = 38.6 

^ (Apf)-mi = -539.93 


Figure 1.- Planform for estimating effect of lee-side pressure on orbiter 
pitching moment. Length in meters; area in m^; planforms indicated 
by numbers. 
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Distance, m 


Figure 2.- Streamline pressure distribution for orbiter configuration starting on 
center line at 9.14 ra from nose, a = 20°; u^ = 6096 m/sec; Altitude = 67 056 m; 
P2 = (3.9 kPa X 10“2 atm); p^ = (0.095 kPa x 10“^ atm). 
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Figure 3.- Isentropic exponent in expansion to leeward side. 
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Figure 4.- Variation of pressure coefficient with flow deflection angle 
Velocity, 6096 m/sec; altitude, 60 964 m. 




Figure 5.- Oblique shock isentropic exponent for 


OrViiter altitude 
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Figure 6.- Isentropic exponent for air. Altitude, 60 960 m, velocity, 6705.6 m/seo 




Figure Real-gas effect on flap effectiveness. 
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Figure 8.- Real-gas effect on pitching- moment coefficient. 
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Figure 10.- Hypersonic viscous interaction parameter behind wedge shock as 
function of Mach number behind wedge shock for shuttle trajectory and 
Langley facilities. 
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Figure 11.- Reynolds number as function of stream Mach number for shuttle 
trajectory and Langley hypersonic facilities. 
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12.- Reynolds number behind wedge shock as function of Mach number behind wedge shock 
for shuttle trajectory (a = 42*^ position) and Langley hypersonic facilities. 
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Figure 15.- Viscous pitching-moment coefficient as function of viscous interaction 

parameter based on free-stream conditions. 




Figure 16.- Viscous drag coefficient as function of viscous interaction 
parameter based on local conditions at x/i^ = O.'i. 



Figure 17.- Viscous pitching moment as function of viscous interaction 
parameter based on local conditions at x/i^ = 0.5. 














Figure 22.- Inviscid drag coefficient on 40^ half-angle cone (obtained with partially coupled 

inviscid-viscid real-gas solution). 
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Figure 24.- Inviscid correlation of pressure 
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Figure 26.- Correlation of total drag coefficient (nondimensionalized by local dynamic pressure) 
with local Mach number. Real-gas effects correlate best in terms of with coefficients 

nondimensionalized by q-^ . 



NASA-Langley, 1977 L-1080? 


I 





Altitude, 

km 

76.2 

70.1 

64.0 

64.0 


Air 

^2 

He 

CF. 


1 


2 



4 


5 


Figure 27.- Correlation of total pitching-moment coefficient (nondimensionalized by local dynamic 

pressure) with local Mach number. 
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